(* Content-type: application/mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 6.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       145,          7]
NotebookDataLength[     71842,       1950]
NotebookOptionsPosition[     66195,       1761]
NotebookOutlinePosition[     66538,       1776]
CellTagsIndexPosition[     66495,       1773]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[{
 RowBox[{
  RowBox[{"Off", "[", 
   RowBox[{"ClebschGordan", "::", "\"\<tri\>\""}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Off", "[", 
   RowBox[{"ClebschGordan", "::", "\"\<phy\>\""}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Off", "[", 
   RowBox[{"SixJSymbol", "::", "\"\<tri\>\""}], "]"}], ";"}]}], "Input",
 CellChangeTimes->{{3.456663504453125*^9, 3.45666350465625*^9}, {
  3.5407344053453703`*^9, 3.540734425171504*^9}}],

Cell[CellGroupData[{

Cell[TextData[{
 "\n",
 StyleBox["Two-photon line shape calculation for K-39 ", "Title"],
 StyleBox["D-2 line from 4S to 4D5/2", "Title",
  FontFamily->"Arial",
  FontVariations->{"Underline"->True}]
}], "Input",
 CellChangeTimes->{{3.551727958766342*^9, 3.5517280116923695`*^9}, {
   3.5517283803094535`*^9, 3.55172838042146*^9}, {3.5521680002246284`*^9, 
   3.552168000672654*^9}, 3.5527533840615234`*^9, {3.5527566752841797`*^9, 
   3.5527566754521484`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"-", "39"}], " ", "D"}], "-", 
  RowBox[{
  "calculation", " ", "for", " ", "K", " ", "line", " ", "photon", " ", 
   "shape"}], "-", 
  RowBox[{"16", " ", "D5", " ", "from", " ", "line", " ", "S", " ", "to"}], 
  "+", "Two"}]], "Output",
 CellChangeTimes->{3.5521685414453125`*^9, 3.552753832883789*^9, 
  3.5527572794716797`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Ouput Location ", "Subsubtitle",
 CellChangeTimes->{{3.5517158329247837`*^9, 3.551715854388012*^9}, 
   3.551716247196479*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{
   "oFile", "=", 
    "\"\<C:\\\\Users\\\\work\\\\Data\\\\output\\\\K39-D2-4D52.txt\>\""}], 
   ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{"output", " ", "file"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.5407383642928085`*^9, 3.5407383651238565`*^9}, 
   3.5407384082223215`*^9, {3.5407384849037075`*^9, 3.540738535987629*^9}, {
   3.54073857516387*^9, 3.5407386491461015`*^9}, {3.5407388655584793`*^9, 
   3.5407388798112946`*^9}, {3.5407453919391146`*^9, 3.5407453981011543`*^9}, 
   3.540745430736563*^9, {3.5409092806939077`*^9, 3.540909299115961*^9}, {
   3.5409139555021634`*^9, 3.5409139572030754`*^9}, {3.540916902879408*^9, 
   3.5409169748426695`*^9}, {3.5409174063863525`*^9, 
   3.5409174253559527`*^9}, {3.5413429497201304`*^9, 
   3.5413429616052127`*^9}, {3.5413431664128094`*^9, 3.541343179302828*^9}, {
   3.542556628837941*^9, 3.542556640226014*^9}, {3.54255670671364*^9, 
   3.5425567165729036`*^9}, {3.5425569015744896`*^9, 3.542556922993427*^9}, {
   3.5425654644126453`*^9, 3.542565470200283*^9}, {3.542570842449442*^9, 
   3.5425708492199287`*^9}, {3.542852782822707*^9, 3.5428527850223217`*^9}, {
   3.542854309978897*^9, 3.542854310868102*^9}, {3.54299342661618*^9, 
   3.5429934460288*^9}, {3.5429966016528063`*^9, 3.5429966440456142`*^9}, 
   3.5429966817575817`*^9, {3.5486039749375*^9, 3.5486040204375*^9}, {
   3.548604193265625*^9, 3.548604193390625*^9}, {3.549060869234375*^9, 
   3.54906087171875*^9}, {3.5492471925820312`*^9, 3.5492471933359375`*^9}, {
   3.5500073476933594`*^9, 3.5500073526621094`*^9}, {3.5512341310498047`*^9, 
   3.5512341453935547`*^9}, {3.5517158046171646`*^9, 3.551715826351408*^9}, 
   3.551716178698561*^9, {3.5517175755924587`*^9, 3.551717580352731*^9}, {
   3.5517284278251715`*^9, 3.551728427945178*^9}, {3.5521680129523563`*^9, 
   3.552168013395382*^9}, {3.552175399694336*^9, 3.552175401084961*^9}, {
   3.5527533777617188`*^9, 3.5527533778847656`*^9}, {3.552756679625*^9, 
   3.5527566798066406`*^9}}],

Cell[TextData[{
 StyleBox["Uncomment the code below to use ",
  FontFamily->"Times New Roman",
  FontWeight->"Bold",
  FontColor->RGBColor[1., 0., 0.]],
 StyleBox["experimental",
  FontFamily->"Times New Roman",
  FontWeight->"Bold",
  FontVariations->{"Underline"->True},
  FontColor->RGBColor[1., 0., 0.]],
 StyleBox[" reptition rate and offset frequency",
  FontFamily->"Times New Roman",
  FontWeight->"Bold",
  FontColor->RGBColor[1., 0., 0.]]
}], "Text",
 CellChangeTimes->{{3.5517161881951046`*^9, 3.551716263157392*^9}, {
  3.5517170355175686`*^9, 3.551717038884761*^9}, {3.5517277525975504`*^9, 
  3.5517277846923857`*^9}},
 FormatType->"TextForm"],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
   RowBox[{
   "pathname", "=", "\"\<C:\\\\Users\\\\work\\\\Data\\\\2012-05-15\\\\\>\""}],
    ";", "\[IndentingNewLine]", 
   RowBox[{"filename", "=", "\"\<Rb-D2-line.dat\>\""}], ";", " ", 
   RowBox[{"(*", "datafile", "*)"}], "\[IndentingNewLine]", 
   RowBox[{"ExpSpecFile", "=", 
    RowBox[{"pathname", " ", "<>", " ", "filename"}]}], ";", 
   "\[IndentingNewLine]", 
   RowBox[{"dataTab", "=", " ", 
    RowBox[{"Import", "[", 
     RowBox[{"ExpSpecFile", ",", "\"\<Table\>\""}], "]"}]}], ";", 
   "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{"Form", " ", "repetition", " ", "rate", " ", 
     RowBox[{"list", ".", " ", "Extract"}], " ", "from", " ", "data", " ", 
     "file"}], "*)"}], "\[IndentingNewLine]", 
   RowBox[{"fOffsetList", "=", 
    RowBox[{"{", "}"}]}], ";", "\[IndentingNewLine]", 
   RowBox[{"fRepList", "=", 
    RowBox[{"{", "}"}]}], ";", "\[IndentingNewLine]", 
   RowBox[{"n", "=", "1"}], ";", "\[IndentingNewLine]", 
   RowBox[{"While", "[", 
    RowBox[{
     RowBox[{"n", "\[LessEqual]", " ", 
      RowBox[{"Length", "[", "dataTab", "]"}]}], ",", "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"fRepList", "=", 
        RowBox[{"Append", "[", 
         RowBox[{"fRepList", ",", 
          RowBox[{
           RowBox[{"dataTab", "[", 
            RowBox[{"[", "n", "]"}], "]"}], "[", 
           RowBox[{"[", "2", "]"}], "]"}]}], "]"}]}], ";", " ", 
       RowBox[{"(*", 
        RowBox[{"\"\<2\>\"", " ", "=", " ", 
         RowBox[{
         "2", "nd", " ", "column", " ", "of", " ", "the", " ", "datafile", 
          " ", "contains", " ", "the", " ", "rep", " ", "rate"}]}], "*)"}], 
       "\[IndentingNewLine]", 
       RowBox[{"fOffsetList", "=", 
        RowBox[{"Append", "[", 
         RowBox[{"fOffsetList", ",", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "2"}], "*", 
            RowBox[{
             RowBox[{"dataTab", "[", 
              RowBox[{"[", "n", "]"}], "]"}], "[", 
             RowBox[{"[", "1", "]"}], "]"}]}], ")"}]}], "]"}]}], ";", " ", 
       RowBox[{"(*", " ", 
        RowBox[{
         RowBox[{"+", 
          RowBox[{"/", 
           RowBox[{"-", " ", "depends"}]}]}], " ", "on", " ", "locking", " ", 
         "electronics"}], "*)"}], "\[IndentingNewLine]", 
       RowBox[{"n", "++"}], ";"}], "}"}]}], "]"}], ";", "\[IndentingNewLine]", 
   RowBox[{"Clear", "[", "n", "]"}], ";", "\[IndentingNewLine]", 
   RowBox[{"fRepMin", "=", 
    RowBox[{"Min", "[", "fRepList", "]"}]}], ";", "\[IndentingNewLine]", 
   RowBox[{"fRepMax", "=", 
    RowBox[{"Max", "[", "fRepList", "]"}]}], ";", "\[IndentingNewLine]", 
   RowBox[{"fRepAvg", "=", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"fRepMin", "+", "fRepMax"}], ")"}], "/", "2."}]}], ";"}], 
  "*)"}]], "Input",
 CellChangeTimes->{{3.5407383642928085`*^9, 3.5407383651238565`*^9}, 
   3.5407384082223215`*^9, {3.5407384849037075`*^9, 3.540738535987629*^9}, {
   3.54073857516387*^9, 3.5407386491461015`*^9}, {3.5407388655584793`*^9, 
   3.5407388798112946`*^9}, {3.5407453919391146`*^9, 3.5407453981011543`*^9}, 
   3.540745430736563*^9, {3.5409092806939077`*^9, 3.540909299115961*^9}, {
   3.5409139555021634`*^9, 3.5409139572030754`*^9}, {3.540916902879408*^9, 
   3.5409169748426695`*^9}, {3.5409174063863525`*^9, 
   3.5409174253559527`*^9}, {3.5413429497201304`*^9, 
   3.5413429616052127`*^9}, {3.5413431664128094`*^9, 3.541343179302828*^9}, {
   3.542556628837941*^9, 3.542556640226014*^9}, {3.54255670671364*^9, 
   3.5425567165729036`*^9}, {3.5425569015744896`*^9, 3.542556922993427*^9}, {
   3.5425654644126453`*^9, 3.542565470200283*^9}, {3.542570842449442*^9, 
   3.5425708492199287`*^9}, {3.542852782822707*^9, 3.5428527850223217`*^9}, {
   3.542854309978897*^9, 3.542854310868102*^9}, {3.54299342661618*^9, 
   3.5429934460288*^9}, {3.5429966016528063`*^9, 3.5429966440456142`*^9}, 
   3.5429966817575817`*^9, {3.5486039749375*^9, 3.5486040204375*^9}, {
   3.548604193265625*^9, 3.548604193390625*^9}, {3.549060869234375*^9, 
   3.54906087171875*^9}, {3.5492471925820312`*^9, 3.5492471933359375`*^9}, {
   3.5500073476933594`*^9, 3.5500073526621094`*^9}, {3.5512341310498047`*^9, 
   3.5512341453935547`*^9}, {3.5517158046171646`*^9, 3.551715826351408*^9}, 
   3.551716178698561*^9, {3.55171684287755*^9, 3.5517168751273947`*^9}, {
   3.5517169260393066`*^9, 3.5517169994325047`*^9}, {3.551727761577064*^9, 
   3.551727761648068*^9}}],

Cell[TextData[{
 StyleBox["Uncomment the code below to use ",
  FontFamily->"Times New Roman",
  FontWeight->"Bold",
  FontColor->RGBColor[1., 0., 0.]],
 StyleBox["theoretical",
  FontFamily->"Times New Roman",
  FontWeight->"Bold",
  FontVariations->{"Underline"->True},
  FontColor->RGBColor[1., 0., 0.]],
 StyleBox[" repitition rate and offset frequency",
  FontFamily->"Times New Roman",
  FontWeight->"Bold",
  FontColor->RGBColor[1., 0., 0.]],
 " "
}], "Text",
 CellChangeTimes->{{3.551716164954775*^9, 3.5517161669018865`*^9}, {
  3.55171627657716*^9, 3.551716297164337*^9}, {3.5517278230845823`*^9, 
  3.5517278563484845`*^9}},
 FormatType->"TextForm"],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"f0", "=", 
    RowBox[{
     RowBox[{"-", "2"}], 
     RowBox[{"(", 
      RowBox[{"20.", " ", 
       SuperscriptBox["10", "6"]}], ")"}]}]}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"after", " ", 
     RowBox[{"SHG", ".", " ", "Add"}]}], " ", "+", 
    RowBox[{"/", "-"}]}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"fRepMin", "=", "250000000."}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{"default", " ", "frep"}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"fRepRange", "=", "80."}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{"scanning", " ", "range"}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fRepMax", "=", 
   RowBox[{"fRepMin", "+", "fRepRange"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Delta]f", "=", "0.25"}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{"Hz", " ", "-", " ", 
    RowBox[{"step", " ", "size", " ", "in", " ", "rep", " ", "rate"}]}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fr", "=", "fRepMin"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fRepList", "=", 
   RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"fOffsetList", "=", 
    RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
   "Form", " ", "the", " ", "fRepList", " ", "and", " ", "fOffsetList"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"While", " ", "[", 
   RowBox[{
    RowBox[{"fr", "<", "fRepMax"}], ",", "\[IndentingNewLine]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"fRepList", "=", 
       RowBox[{"Append", "[", 
        RowBox[{"fRepList", ",", "fr"}], "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"fOffsetList", "=", 
       RowBox[{"Append", "[", 
        RowBox[{"fOffsetList", ",", "f0"}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"fr", "+=", "\[Delta]f"}], ";"}], "}"}]}], "]"}], 
  "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{"Average", " ", "fRep"}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fRepAvg", "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"fRepMin", "+", "fRepMax"}], ")"}], "/", "2."}]}], 
  ";"}]}], "Input",
 CellChangeTimes->{{3.5517161466957307`*^9, 3.5517161587894225`*^9}, {
   3.5517163024016366`*^9, 3.5517163047197695`*^9}, {3.5517163366095934`*^9, 
   3.5517163419118967`*^9}, {3.55171653062969*^9, 3.5517167452519665`*^9}, {
   3.5517167950658154`*^9, 3.551716838278287*^9}, 3.5517170583998775`*^9, {
   3.5517286701030283`*^9, 3.551728670279039*^9}, 3.5521680173266068`*^9, {
   3.5521685349296875`*^9, 3.5521685350234375`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Physical Constants", "Subtitle",
 FontColor->RGBColor[0, 0.356863, 0]],

Cell["Speed of Light", "Subsubtitle"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"c", "=", "29979245800."}], ";"}], 
  RowBox[{"(*", 
   RowBox[{"cm", "/", "sec"}], "*)"}]}]], "Input",
 CellChangeTimes->{3.551728122119686*^9}],

Cell[CellGroupData[{

Cell["Planck's Constant", "Subsubtitle"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"hBar", "=", 
    RowBox[{"1.054", " ", 
     SuperscriptBox["10", 
      RowBox[{"-", "27"}]]}]}], ";"}], 
  RowBox[{"(*", 
   RowBox[{"ergs", " ", "sec"}], "*)"}]}]], "Input"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Boltzmann Constant", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"kB", "=", 
    RowBox[{"1.380658", " ", 
     SuperscriptBox["10", 
      RowBox[{"-", "16"}]]}]}], ";"}], 
  RowBox[{"(*", 
   RowBox[{"erg", "/", "K"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.4566596001875*^9, 3.456659625953125*^9}, {
  3.45665968075*^9, 3.456659686359375*^9}}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Atomic Parameters", "Subtitle",
 FontColor->RGBColor[0, 0.356863, 0]],

Cell["General Parameters", "Subtitle",
 CellChangeTimes->{{3.45666146840625*^9, 3.45666147296875*^9}}],

Cell["Nuclear Spin", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"i", "=", 
   RowBox[{"N", "[", 
    FractionBox["3", "2"], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.4730885615*^9, 3.473088561953125*^9}, 3.53632801824225*^9, {
   3.5425567297081876`*^9, 3.5425567322042036`*^9}, {3.5512806388623047`*^9, 
   3.5512806390654297`*^9}}],

Cell[CellGroupData[{

Cell["Transit Linewidth", "Subsubtitle",
 CellChangeTimes->{{3.45666389359375*^9, 3.456663897328125*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"\[Gamma]T", "=", 
   FractionBox["1", 
    RowBox[{"40.2", " ", 
     SuperscriptBox["10", 
      RowBox[{"-", "9"}]]}]]}], ";"}]], "Input",
 CellChangeTimes->{{3.456663726421875*^9, 3.4566637388125*^9}, {
   3.45666380203125*^9, 3.456663813921875*^9}, 3.4566922269375*^9, 
   3.5366835360235*^9, 3.538523997111776*^9, {3.5492471882460938`*^9, 
   3.549247189754883*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Polarization", "Subsubtitle",
 CellChangeTimes->{{3.456663909203125*^9, 3.456663910671875*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"q1", "=", 
   RowBox[{"{", 
    RowBox[{"0", ",", "1", ",", "0"}], "}"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.456663912765625*^9, 3.456663917453125*^9}, 
   3.456663957921875*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"q2", "=", 
   RowBox[{"{", 
    RowBox[{"0", ",", "1", ",", "0"}], "}"}]}], ";"}]], "Input",
 CellChangeTimes->{3.4566639629375*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Ground State Parameters", "Subtitle",
 CellChangeTimes->{{3.456659852765625*^9, 3.456659860875*^9}}],

Cell["Ground State S", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
   3.45665970596875*^9, 3.456659711984375*^9}, {3.45665974584375*^9, 
   3.456659752203125*^9}, 3.45665980753125*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"SGnd", "=", 
   RowBox[{"N", "[", 
    FractionBox["1", "2"], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.45665975521875*^9, 3.456659757859375*^9}, 3.4566598095625*^9, {
   3.5425567378202395`*^9, 3.5425567406126575`*^9}}],

Cell[CellGroupData[{

Cell["Ground State J", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.45665974584375*^9, 
  3.456659752203125*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"jGnd", "=", 
   RowBox[{"N", "[", 
    FractionBox["1", "2"], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
  3.45665975521875*^9, 3.456659757859375*^9}, {3.5425567479135046`*^9, 
  3.542556749863517*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Ground State F Values", "Subsubtitle",
 CellChangeTimes->{{3.45666290321875*^9, 3.456662908015625*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"fMinGnd", "=", 
   RowBox[{"Abs", "[", 
    RowBox[{"jGnd", "-", "i"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fMaxGnd", "=", 
   RowBox[{"jGnd", "+", "i"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.456662861609375*^9, 3.456662890515625*^9}, {
   3.54906118934375*^9, 3.549061201671875*^9}, {3.5490736629375*^9, 
   3.549073666609375*^9}, {3.5492478572304688`*^9, 3.549247859171875*^9}, {
   3.550007296041992*^9, 3.5500073055283203`*^9}, {3.5500073762246094`*^9, 
   3.5500073772871094`*^9}, 3.5512806623154297`*^9, {3.5512807609091797`*^9, 
   3.5512807639716797`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Ground State Hyperfine A", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
  3.45665976646875*^9}, {3.456659909703125*^9, 3.45665991996875*^9}, {
  3.45666013296875*^9, 3.456660152140625*^9}, {3.456660398109375*^9, 
  3.45666040421875*^9}, {3.456660756828125*^9, 3.45666075759375*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"aGnd", "=", 
    RowBox[{"SetPrecision", "[", 
     RowBox[{"230859860100.", ",", "20"}], "]"}]}], " ", ";"}], 
  RowBox[{"(*", " ", 
   RowBox[{"230.8598601", " ", "MHz"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.45665976928125*^9, 3.45665977578125*^9}, {3.45665992590625*^9, 
   3.45665992978125*^9}, {3.456659964359375*^9, 3.45665998015625*^9}, {
   3.456660078296875*^9, 3.456660092390625*^9}, {3.456660306609375*^9, 
   3.45666036596875*^9}, {3.456660407359375*^9, 3.45666043703125*^9}, {
   3.4566605581875*^9, 3.4566605719375*^9}, {3.456660680015625*^9, 
   3.456660705609375*^9}, {3.45666076115625*^9, 3.4566607615625*^9}, {
   3.47308857578125*^9, 3.473088584546875*^9}, {3.47309640846875*^9, 
   3.473096424140625*^9}, {3.536335670945375*^9, 3.53633567542975*^9}, {
   3.5512809188779297`*^9, 3.5512809360341797`*^9}, {3.5512810990966797`*^9, 
   3.5512811179873047`*^9}, 3.5512811900185547`*^9, {3.551715895359355*^9, 
   3.5517159071990323`*^9}, {3.551717675681184*^9, 3.5517176816565256`*^9}, 
   3.551728139904703*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Ground State Hyperfine B", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
   3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
   3.45665976646875*^9}, {3.456659909703125*^9, 3.45665991996875*^9}, {
   3.45666013296875*^9, 3.456660152140625*^9}, {3.456660398109375*^9, 
   3.45666040421875*^9}, 3.4566605775*^9, {3.45666076390625*^9, 
   3.45666076546875*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"bGnd", "=", "0."}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.45665976928125*^9, 3.45665977578125*^9}, {3.45665992590625*^9, 
   3.45665992978125*^9}, {3.456659964359375*^9, 3.45665998015625*^9}, {
   3.456660078296875*^9, 3.456660092390625*^9}, {3.456660306609375*^9, 
   3.45666036596875*^9}, {3.456660407359375*^9, 3.45666043703125*^9}, {
   3.4566605581875*^9, 3.456660602203125*^9}, {3.45666063840625*^9, 
   3.45666063971875*^9}, {3.45666076803125*^9, 3.456660773828125*^9}, 
   3.5512811922998047`*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[TextData[{
 "Intermediate State Parameters ",
 Cell[BoxData[
  FormBox[
   RowBox[{"4", 
    SubscriptBox["P", 
     RowBox[{"3", "/", "2"}]]}], TraditionalForm]]]
}], "Subtitle",
 CellChangeTimes->{{3.456659852765625*^9, 3.456659887765625*^9}, {
  3.5512809799404297`*^9, 3.5512809843154297`*^9}, {3.5527533931523438`*^9, 
  3.5527533954726562`*^9}}],

Cell["Intermediate State S", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.45665974584375*^9, 
  3.456659752203125*^9}, {3.45665980753125*^9, 3.45665981840625*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"SInt", "=", 
   RowBox[{"N", "[", 
    FractionBox["1", "2"], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
  3.45665975521875*^9, 3.456659757859375*^9}, {3.4566598095625*^9, 
  3.45665982178125*^9}, {3.542556758911575*^9, 3.5425567611423893`*^9}}],

Cell[CellGroupData[{

Cell["Intermediate State J", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
  3.45665976646875*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"jInt", "=", 
   RowBox[{"N", "[", 
    FractionBox["3", "2"], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.45665976928125*^9, 3.45665977578125*^9}, {3.473088588953125*^9, 
   3.47308858946875*^9}, 3.536333582476625*^9, {3.5425567651828156`*^9, 
   3.5425567679908333`*^9}, {3.5512809908779297`*^9, 
   3.5512809910029297`*^9}, {3.5527533995458984`*^9, 3.5527533996796875`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Intermediate State F Values", "Subsubtitle",
 CellChangeTimes->{{3.45666290321875*^9, 3.456662923203125*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"fMinInt", "=", 
   RowBox[{"Abs", "[", 
    RowBox[{"jInt", "-", "i"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fMaxInt", "=", 
   RowBox[{"jInt", "+", "i"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.456662861609375*^9, 3.456662890515625*^9}, {
  3.456662927125*^9, 3.456662936640625*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Intermediate State Linewidth", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
  3.45665976646875*^9}, {3.456659909703125*^9, 3.45665991996875*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"\[Gamma]Int", "=", 
    RowBox[{"N", "[", 
     FractionBox["1", 
      RowBox[{"26.37", " ", 
       SuperscriptBox["10", 
        RowBox[{"-", "9"}]]}]], "]"}]}], ";"}], 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"Seconds", " ", "From", " ", "Budker"}], ",", "Kimball", ",", " ",
     "DeMille"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
  3.45665976928125*^9, 3.45665977578125*^9}, {3.45665992590625*^9, 
  3.45665992978125*^9}, {3.456659964359375*^9, 3.45665998015625*^9}, {
  3.456660078296875*^9, 3.456660092390625*^9}, {3.53668310449225*^9, 
  3.536683106007875*^9}, {3.5512810147529297`*^9, 3.5512810175185547`*^9}, {
  3.5512811805185547`*^9, 3.5512811853779297`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Intermediate State Center of Gravity Energy", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
  3.45665976646875*^9}, {3.456659909703125*^9, 3.45665991996875*^9}, {
  3.45666013296875*^9, 3.456660152140625*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"\[Nu]Int", "=", 
   RowBox[{"SetPrecision", "[", 
    RowBox[{
     RowBox[{"391.01617083", " ", 
      SuperscriptBox["10", "12"]}], ",", "20"}], "]"}]}], " ", 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"Hz", " ", "from", " ", "Steck"}], "*)"}]}]}], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.45665976928125*^9, 3.45665977578125*^9}, {3.45665992590625*^9, 
   3.45665992978125*^9}, {3.456659964359375*^9, 3.45665998015625*^9}, {
   3.456660078296875*^9, 3.456660092390625*^9}, {3.456660306609375*^9, 
   3.45666036596875*^9}, {3.473088596515625*^9, 3.473088616875*^9}, {
   3.47308883840625*^9, 3.473088843484375*^9}, {3.473096433640625*^9, 
   3.473096439859375*^9}, {3.53633478836725*^9, 3.53633479311725*^9}, 
   3.536607319211*^9, {3.536607447226625*^9, 3.536607498007875*^9}, {
   3.53666646042975*^9, 3.5366664651485*^9}, {3.537893167558856*^9, 
   3.5378931893520937`*^9}, 3.5378943549667425`*^9, 3.5378974743522215`*^9, {
   3.5378985672711487`*^9, 3.53789856825395*^9}, {3.5378998431661897`*^9, 
   3.537899869374236*^9}, 3.5379017571978636`*^9, 3.538606320562795*^9, {
   3.5512810393310547`*^9, 3.5512810634560547`*^9}, 3.5512811714248047`*^9, {
   3.552753582833008*^9, 3.552753602392578*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Intermediate State Hyperfine A", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
  3.45665976646875*^9}, {3.456659909703125*^9, 3.45665991996875*^9}, {
  3.45666013296875*^9, 3.456660152140625*^9}, {3.456660398109375*^9, 
  3.45666040421875*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"aInt", "=", 
    RowBox[{"SetPrecision", "[", 
     RowBox[{
      RowBox[{"6.933", " ", 
       SuperscriptBox["10", "6"]}], ",", "20"}], "]"}]}], " ", ";"}], 
  RowBox[{"(*", 
   RowBox[{"Hz", " ", "From", " ", "Steck"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
  3.45665976928125*^9, 3.45665977578125*^9}, {3.45665992590625*^9, 
  3.45665992978125*^9}, {3.456659964359375*^9, 3.45665998015625*^9}, {
  3.456660078296875*^9, 3.456660092390625*^9}, {3.456660306609375*^9, 
  3.45666036596875*^9}, {3.456660407359375*^9, 3.45666043703125*^9}, {
  3.4566605581875*^9, 3.4566605719375*^9}, {3.473088627078125*^9, 
  3.4730886284375*^9}, {3.473096443796875*^9, 3.473096448328125*^9}, {
  3.5363357313985*^9, 3.536335734726625*^9}, {3.536607079039125*^9, 
  3.53660709955475*^9}, {3.536666472757875*^9, 3.536666476257875*^9}, {
  3.536683146882875*^9, 3.5366831523985*^9}, {3.5372022874831247`*^9, 
  3.5372022919049997`*^9}, {3.5512811388466797`*^9, 3.5512811690185547`*^9}, {
  3.5517160806559534`*^9, 3.551716087167326*^9}, {3.5527536052421875`*^9, 
  3.5527536164472656`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Intermediate State Hyperfine B", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
   3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
   3.45665976646875*^9}, {3.456659909703125*^9, 3.45665991996875*^9}, {
   3.45666013296875*^9, 3.456660152140625*^9}, {3.456660398109375*^9, 
   3.45666040421875*^9}, 3.4566605775*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"bInt", "=", 
    RowBox[{"SetPrecision", "[", 
     RowBox[{
      RowBox[{"2.786", " ", 
       SuperscriptBox["10", "6"]}], ",", "20"}], "]"}]}], " ", ";"}], 
  RowBox[{"(*", 
   RowBox[{"Hz", " ", "From", " ", "Steck"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.45665976928125*^9, 3.45665977578125*^9}, {3.45665992590625*^9, 
   3.45665992978125*^9}, {3.456659964359375*^9, 3.45665998015625*^9}, {
   3.456660078296875*^9, 3.456660092390625*^9}, {3.456660306609375*^9, 
   3.45666036596875*^9}, {3.456660407359375*^9, 3.45666043703125*^9}, {
   3.4566605581875*^9, 3.456660602203125*^9}, {3.45666063840625*^9, 
   3.45666063971875*^9}, 3.473088631296875*^9, {3.473096452546875*^9, 
   3.473096456171875*^9}, {3.536335742211*^9, 3.536335745961*^9}, {
   3.536607112711*^9, 3.536607116664125*^9}, 3.53666648311725*^9, {
   3.5366831642735*^9, 3.53668317030475*^9}, {3.5372023036862497`*^9, 
   3.5372023085924997`*^9}, {3.5512811663310547`*^9, 
   3.5512811667529297`*^9}, {3.55275362209375*^9, 3.5527536285117188`*^9}}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[TextData[{
 "Final State Parameters ",
 Cell[BoxData[
  FormBox[
   RowBox[{"4", 
    SubscriptBox["D", 
     RowBox[{"5", "/", "2"}]]}], TraditionalForm]]]
}], "Subtitle",
 CellChangeTimes->{{3.456659852765625*^9, 3.456659887765625*^9}, {
  3.456660101109375*^9, 3.4566601018125*^9}, {3.5512819261904297`*^9, 
  3.5512819303154297`*^9}, {3.551717330572445*^9, 3.5517173308154583`*^9}, {
  3.5517283754521756`*^9, 3.5517283755511813`*^9}, {3.552168156124545*^9, 
  3.552168156534569*^9}, {3.552756769125*^9, 3.5527567692197266`*^9}}],

Cell["Final State S", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.45665974584375*^9, 
  3.456659752203125*^9}, {3.45665980753125*^9, 3.456659829375*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"SFin", "=", 
   RowBox[{"N", "[", 
    FractionBox["1", "2"], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.45665975521875*^9, 3.456659757859375*^9}, {3.4566598095625*^9, 
   3.45665983434375*^9}, {3.536335431601625*^9, 3.536335431726625*^9}, 
   3.53634248861725*^9, {3.5425567824365253`*^9, 3.5425567858841476`*^9}}],

Cell[CellGroupData[{

Cell["Final State J", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.456659783546875*^9, 
  3.456659787671875*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"jFin", "=", 
   RowBox[{"N", "[", 
    FractionBox["5", "2"], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.4566597779375*^9, 3.456659780921875*^9}, {3.456661545015625*^9, 
   3.456661545171875*^9}, 3.4569295484212503`*^9, {3.47308863565625*^9, 
   3.473088640640625*^9}, {3.536335439086*^9, 3.536335439226625*^9}, {
   3.542556787927761*^9, 3.542556790751379*^9}, {3.5512819388935547`*^9, 
   3.5512819389873047`*^9}, {3.5517284386257887`*^9, 3.551728438744796*^9}, {
   3.5521681719094477`*^9, 3.5521681723134713`*^9}, {3.552756771123047*^9, 
   3.552756771241211*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Final State F Values", "Subsubtitle",
 CellChangeTimes->{{3.45666290321875*^9, 3.45666294725*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"fMinFin", "=", 
   RowBox[{"Abs", "[", 
    RowBox[{"jFin", "-", "i"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fMaxFin", "=", 
   RowBox[{"jFin", "+", "i"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.456662861609375*^9, 3.456662890515625*^9}, {
  3.456662927125*^9, 3.456662957578125*^9}, {3.54906121109375*^9, 
  3.54906122009375*^9}, {3.5490736735625*^9, 3.5490736805625*^9}, {
  3.549247847415039*^9, 3.549247849538086*^9}, {3.5500073128964844`*^9, 
  3.5500073230839844`*^9}, {3.5500073833496094`*^9, 3.5500073842246094`*^9}, {
  3.5512819426591797`*^9, 3.5512819468466797`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"\[Nu]Fin", "=", 
   RowBox[{"SetPrecision", "[", 
    RowBox[{
     RowBox[{"27397.077", " ", "c"}], ",", "20"}], "]"}]}], " ", 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"aFin", "=", 
   RowBox[{"SetPrecision", "[", 
    RowBox[{
     RowBox[{"0.120", " ", 
      SuperscriptBox["10", "6"]}], ",", "20"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"bFin", "=", "0.0"}], ";"}]}], "Input",
 CellChangeTimes->{{3.5512838488154297`*^9, 3.5512838633310547`*^9}, {
   3.551717423621767*^9, 3.551717433431328*^9}, 3.5517275821998043`*^9, 
   3.5517281530804567`*^9, {3.5517283163917975`*^9, 3.551728322063122*^9}, {
   3.55172845647081*^9, 3.5517284682794847`*^9}, {3.552168488942581*^9, 
   3.552168492693796*^9}, {3.5527572250371094`*^9, 3.5527572423964844`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Final State Linewidth", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
  3.45665976646875*^9}, {3.456659909703125*^9, 3.45665991996875*^9}, {
  3.456659988796875*^9, 3.45665998953125*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"\[Gamma]Fin", "=", 
    FractionBox["1", 
     RowBox[{"291.18", " ", 
      SuperscriptBox["10", 
       RowBox[{"-", "9"}]]}]]}], ";"}], 
  RowBox[{"(*", 
   RowBox[{"Seconds", " ", "From", " ", "Heavens"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
  3.45665976928125*^9, 3.45665977578125*^9}, {3.45665992590625*^9, 
  3.45665992978125*^9}, {3.456659964359375*^9, 3.45665998015625*^9}, {
  3.4566600681875*^9, 3.45666007421875*^9}, {3.4566616303125*^9, 
  3.456661631765625*^9}, {3.4569295617493753`*^9, 3.4569295620150003`*^9}, {
  3.45951917796875*^9, 3.45951917865625*^9}, {3.536683236757875*^9, 
  3.536683244726625*^9}, {3.538582219701305*^9, 3.5385822378273416`*^9}, {
  3.5512838421123047`*^9, 3.5512838442841797`*^9}, {3.5517274444239235`*^9, 
  3.5517274454019794`*^9}, {3.551727533652027*^9, 3.5517275372902355`*^9}, {
  3.5521683695107503`*^9, 3.5521683771941895`*^9}, {3.552757252254883*^9, 
  3.552757258881836*^9}}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Energies, Matrix Elements and Transition Probablities", "Subtitle",
 CellChangeTimes->{{3.456681533796875*^9, 3.456681554234375*^9}}],

Cell["Formula for Calculating Energies with HFS", "Subsubtitle",
 CellChangeTimes->{{3.45666156509375*^9, 3.45666157378125*^9}, {
  3.456662570765625*^9, 3.4566625805625*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"\[Kappa]", "[", 
    RowBox[{"f_", ",", "j_"}], "]"}], "=", 
   RowBox[{
    RowBox[{"f", 
     RowBox[{"(", 
      RowBox[{"f", "+", "1"}], ")"}]}], "-", 
    RowBox[{"j", 
     RowBox[{"(", 
      RowBox[{"j", "+", "1"}], ")"}]}], "-", 
    RowBox[{"i", 
     RowBox[{"(", 
      RowBox[{"i", "+", "1"}], ")"}]}]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.45666169846875*^9, 3.456661713*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"\[Delta]E", "[", 
    RowBox[{"a_", ",", "b_", ",", "f_", ",", "j_"}], "]"}], ":=", 
   RowBox[{"If", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"j", "\[NotEqual]", " ", 
       FractionBox["1", "2"]}], "&&", 
      RowBox[{"i", "\[NotEqual]", 
       FractionBox["1", "2"]}]}], ",", 
     RowBox[{
      RowBox[{
       FractionBox["1", "2"], "a", " ", 
       RowBox[{"\[Kappa]", "[", 
        RowBox[{"f", ",", "j"}], "]"}]}], "+", 
      RowBox[{"b", " ", 
       FractionBox[
        RowBox[{
         RowBox[{
          FractionBox["3", "4"], 
          SuperscriptBox[
           RowBox[{"\[Kappa]", "[", 
            RowBox[{"f", ",", "j"}], "]"}], "2"]}], "+", 
         RowBox[{
          FractionBox["3", "4"], 
          RowBox[{"\[Kappa]", "[", 
           RowBox[{"f", ",", "j"}], "]"}]}], "-", 
         RowBox[{"i", 
          RowBox[{"(", 
           RowBox[{"i", "+", "1"}], ")"}], "j", 
          RowBox[{"(", 
           RowBox[{"j", "+", "1"}], ")"}]}]}], 
        RowBox[{"2", " ", "i", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"2", "i"}], "-", "1"}], ")"}], "j", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"2", "j"}], "-", "1"}], ")"}]}]]}]}], ",", 
     RowBox[{
      FractionBox["1", "2"], "a", " ", 
      RowBox[{"\[Kappa]", "[", 
       RowBox[{"f", ",", "j"}], "]"}]}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.456661820734375*^9, 3.45666183171875*^9}, {
   3.456662393875*^9, 3.456662397109375*^9}, {3.456683422046875*^9, 
   3.456683431546875*^9}, 3.460223934290625*^9}],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"\[Nu]GndToInt", "[", 
    RowBox[{"f_", ",", "fPrime_"}], "]"}], "=", 
   RowBox[{"\[Nu]Int", "-", 
    RowBox[{"\[Delta]E", "[", 
     RowBox[{"aGnd", ",", "bGnd", ",", "f", ",", "jGnd"}], "]"}], "+", 
    RowBox[{"\[Delta]E", "[", 
     RowBox[{"aInt", ",", "bInt", ",", "fPrime", ",", "jInt"}], "]"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Nu]GndToFin", "[", 
    RowBox[{"f_", ",", "fDPrime_"}], "]"}], "=", 
   RowBox[{"\[Nu]Fin", "-", 
    RowBox[{"\[Delta]E", "[", 
     RowBox[{"aGnd", ",", "bGnd", ",", "f", ",", "jGnd"}], "]"}], "+", 
    RowBox[{"\[Delta]E", "[", 
     RowBox[{"aFin", ",", "bFin", ",", "fDPrime", ",", "jFin"}], "]"}]}]}], 
  ";"}]}], "Input",
 CellChangeTimes->{{3.45686233334375*^9, 3.4568623374375*^9}, 
   3.4572661415462503`*^9, {3.47308885534375*^9, 3.473088860125*^9}, {
   3.4730964853125*^9, 3.473096517875*^9}, {3.473097064765625*^9, 
   3.47309707559375*^9}, {3.536683396664125*^9, 3.53668342911725*^9}, {
   3.5385822740464134`*^9, 3.5385822762455387`*^9}, {3.538582433139513*^9, 
   3.538582435955674*^9}, {3.53858286930646*^9, 3.538583146448312*^9}, {
   3.5399865266491117`*^9, 3.539986527335513*^9}, {3.5425567991442327`*^9, 
   3.5425568137147264`*^9}, {3.5517174795439653`*^9, 
   3.5517175218973875`*^9}, {3.551727944195509*^9, 3.551727944470525*^9}, {
   3.551728158825785*^9, 3.5517281597138357`*^9}, {3.552168184344159*^9, 
   3.5521681859422507`*^9}, 3.552168222350333*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{"\[Nu]GndToFin", "[", 
   RowBox[{"2", ",", "3"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.5517281719065332`*^9, 3.551728192390705*^9}}],

Cell[BoxData["8.211705606594516`*^14"], "Output",
 CellChangeTimes->{{3.551728180674035*^9, 3.5517281927097235`*^9}, 
   3.5521685420859375`*^9, 3.552753833303711*^9, 3.55275728003125*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Reduced Matrix Elements in J Basis", "Subsubtitle"],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"dRedJGndToInt", "=", "1"}], ";"}], 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
     RowBox[{"(", "4.51", ")"}], "1.28", " ", 
     SuperscriptBox["10", "6"]}], " ", 
    RowBox[{"(*", 
     RowBox[{"Hz", "/", 
      RowBox[{"(", 
       RowBox[{"V", "/", "cm"}], ")"}]}], "*)"}], ";"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"dRedJIntToFin", "=", "1"}], ";"}], 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
     FractionBox["1", 
      RowBox[{"e", " ", "a0"}]], 
     RowBox[{"\[Sqrt]", 
      RowBox[{"(", 
       RowBox[{
        FractionBox[
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"2", " ", 
             FractionBox["1", "2"]}], "+", "1"}], " ", ")"}], "3", " ", 
          "hBar", " ", 
          SuperscriptBox["c", "3"]}], 
         RowBox[{"4", " ", 
          SuperscriptBox[
           RowBox[{"(", 
            RowBox[{"2", " ", "\[Pi]", " ", "\[Nu]6P12To8S"}], ")"}], "3"]}]],
         "2.04", " ", 
        SuperscriptBox["10", "6"]}], ")"}]}], "1.28", " ", 
     SuperscriptBox["10", "6"]}], " ", 
    RowBox[{"(*", 
     RowBox[{"Hz", "/", 
      RowBox[{"(", 
       RowBox[{"V", "/", "cm"}], ")"}]}], "*)"}], ";"}], "*)"}]}]}], "Input",
 CellChangeTimes->{{3.4566056364375*^9, 3.4566056744375*^9}, {
  3.456662686453125*^9, 3.456662724453125*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Reduced Matrix Elements in F Basis", "Subsubtitle",
 CellChangeTimes->{3.4566632656875*^9}],

Cell[BoxData[
 RowBox[{"Do", "[", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"Do", "[", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      RowBox[{
       SubscriptBox["dRedGndToInt", 
        RowBox[{"f", ",", "fPrime"}]], "=", 
       RowBox[{"\[Sqrt]", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"2", " ", "fPrime"}], "+", "1"}], ")"}], 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"2", "f"}], "+", "1"}], ")"}], 
          SuperscriptBox[
           RowBox[{"SixJSymbol", "[", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{"jInt", ",", "fPrime", ",", "i"}], "}"}], ",", 
             RowBox[{"{", 
              RowBox[{"f", ",", "jGnd", ",", "1"}], "}"}]}], "]"}], "2"], 
          SuperscriptBox["dRedJGndToInt", "2"]}], ")"}]}]}], ";"}], ",", 
     RowBox[{"{", 
      RowBox[{"fPrime", ",", "fMinInt", ",", "fMaxInt"}], "}"}]}], "]"}], ",",
    "\[IndentingNewLine]", 
   RowBox[{"{", 
    RowBox[{"f", ",", "fMinGnd", ",", "fMaxGnd"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.456662795671875*^9, 3.456662835546875*^9}, {
  3.456663002453125*^9, 3.456663020171875*^9}, {3.5372063527331247`*^9, 
  3.5372063614518747`*^9}}],

Cell[BoxData[
 RowBox[{"Do", "[", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"Do", "[", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      SubscriptBox["dRedIntToFin", 
       RowBox[{"fPrime", ",", "fDPrime"}]], "=", 
      RowBox[{"\[Sqrt]", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"2", " ", "fDPrime"}], "+", "1"}], ")"}], 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"2", "fPrime"}], "+", "1"}], ")"}], 
         SuperscriptBox[
          RowBox[{"SixJSymbol", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{"jFin", ",", "fDPrime", ",", "i"}], "}"}], ",", 
            RowBox[{"{", 
             RowBox[{"fPrime", ",", "jInt", ",", "1"}], "}"}]}], "]"}], "2"], 
         SuperscriptBox["dRedJIntToFin", "2"]}], ")"}]}]}], ",", 
     RowBox[{"{", 
      RowBox[{"fDPrime", ",", "fMinFin", ",", "fMaxFin"}], "}"}]}], "]"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{"fPrime", ",", "fMinInt", ",", "fMaxInt"}], "}"}]}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.456663051265625*^9, 3.456663093375*^9}, {
   3.456663125046875*^9, 3.456663154296875*^9}, 3.456663192609375*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Matrix Elements in F Basis", "Subsubtitle",
 CellChangeTimes->{{3.456663276328125*^9, 3.456663278140625*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Do", "[", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Do", "[", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"Do", "[", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"Do", "[", "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{
           SubscriptBox["dGndToInt", 
            RowBox[{"f", ",", "fPrime", ",", "m", ",", "mPrime"}]], "=", 
           RowBox[{
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{"-", "1"}], ")"}], 
             RowBox[{"fPrime", "-", "mPrime"}]], 
            RowBox[{
             UnderoverscriptBox["\[Sum]", 
              RowBox[{"n", "=", "1"}], "3"], 
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                RowBox[{
                "q1", "\[LeftDoubleBracket]", "n", "\[RightDoubleBracket]"}], 
                RowBox[{"ThreeJSymbol", "[", 
                 RowBox[{
                  RowBox[{"{", 
                   RowBox[{"fPrime", ",", 
                    RowBox[{"-", "mPrime"}]}], "}"}], ",", 
                  RowBox[{"{", 
                   RowBox[{"1", ",", 
                    RowBox[{"n", "-", "2"}]}], "}"}], ",", 
                  RowBox[{"{", 
                   RowBox[{"f", ",", "m"}], "}"}]}], "]"}]}], ")"}], 
              SubscriptBox["dRedGndToInt", 
               RowBox[{"f", ",", "fPrime"}]]}]}]}]}], ",", 
          "\[IndentingNewLine]", 
          RowBox[{"{", 
           RowBox[{"m", ",", 
            RowBox[{"-", "f"}], ",", "f"}], "}"}]}], "]"}], ",", 
        "\[IndentingNewLine]", 
        RowBox[{"{", 
         RowBox[{"mPrime", ",", 
          RowBox[{"-", "fPrime"}], ",", "fPrime"}], "}"}]}], "]"}], ",", 
      "\[IndentingNewLine]", 
      RowBox[{"{", 
       RowBox[{"fPrime", ",", "fMinInt", ",", "fMaxInt"}], "}"}]}], "]"}], 
    ",", "\[IndentingNewLine]", 
    RowBox[{"{", 
     RowBox[{"f", ",", "fMinGnd", ",", "fMaxGnd"}], "}"}]}], "]"}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.45666329196875*^9, 3.4566633891875*^9}, {3.456663993*^9,
   3.45666408128125*^9}, {3.45666411909375*^9, 3.456664121515625*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Do", "[", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Do", "[", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"Do", "[", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"Do", "[", "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{
           SubscriptBox["dIntToFin", 
            RowBox[{
            "fPrime", ",", "fDPrime", ",", "mPrime", ",", "mDPrime"}]], "=", 
           RowBox[{
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{"-", "1"}], ")"}], 
             RowBox[{"fDPrime", "-", "mDPrime"}]], 
            RowBox[{
             UnderoverscriptBox["\[Sum]", 
              RowBox[{"n", "=", "1"}], "3"], 
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                RowBox[{
                "q2", "\[LeftDoubleBracket]", "n", "\[RightDoubleBracket]"}], 
                RowBox[{"ThreeJSymbol", "[", 
                 RowBox[{
                  RowBox[{"{", 
                   RowBox[{"fDPrime", ",", 
                    RowBox[{"-", "mDPrime"}]}], "}"}], ",", 
                  RowBox[{"{", 
                   RowBox[{"1", ",", 
                    RowBox[{"n", "-", "2"}]}], "}"}], ",", 
                  RowBox[{"{", 
                   RowBox[{"fPrime", ",", "mPrime"}], "}"}]}], "]"}]}], ")"}], 
              SubscriptBox["dRedIntToFin", 
               RowBox[{"fPrime", ",", "fDPrime"}]]}]}]}]}], ",", 
          "\[IndentingNewLine]", 
          RowBox[{"{", 
           RowBox[{"mPrime", ",", 
            RowBox[{"-", "fPrime"}], ",", "fPrime"}], "}"}]}], "]"}], ",", 
        "\[IndentingNewLine]", 
        RowBox[{"{", 
         RowBox[{"mDPrime", ",", 
          RowBox[{"-", "fDPrime"}], ",", "fDPrime"}], "}"}]}], "]"}], ",", 
      "\[IndentingNewLine]", 
      RowBox[{"{", 
       RowBox[{"fDPrime", ",", "fMinFin", ",", "fMaxFin"}], "}"}]}], "]"}], 
    ",", "\[IndentingNewLine]", 
    RowBox[{"{", 
     RowBox[{"fPrime", ",", "fMinInt", ",", "fMaxInt"}], "}"}]}], "]"}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.45666329196875*^9, 3.456663489203125*^9}, {
  3.456664105703125*^9, 3.45666415221875*^9}}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Doppler Distribution Parameters", "Subtitle",
 FontColor->RGBColor[0, 0.356863, 0]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"amu", "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"1.66053886", " ", 
      SuperscriptBox["10", 
       RowBox[{"-", "27"}]]}], ")"}], 
    SuperscriptBox["10", "3"]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"mK", "=", 
   RowBox[{"38.96370668", "amu"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"temp", "=", "343"}], ";"}]}], "Input",
 CellChangeTimes->{{3.456606503203125*^9, 3.4566065051875*^9}, {
   3.456664763890625*^9, 3.456664790515625*^9}, {3.473097352640625*^9, 
   3.473097368078125*^9}, {3.536605457211*^9, 3.53660545749225*^9}, 
   3.5366060920235*^9, {3.536606131101625*^9, 3.536606151414125*^9}, 
   3.53666517880475*^9, {3.540734490795257*^9, 3.5407345006498213`*^9}, {
   3.5407345404901*^9, 3.5407345411041346`*^9}, {3.5517172334698906`*^9, 
   3.5517172625375533`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"\[Alpha]", "=", 
   RowBox[{"\[Sqrt]", 
    RowBox[{"(", 
     FractionBox[
      RowBox[{"2", " ", "kB", " ", "temp"}], "mK"], ")"}]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.45666479665625*^9, 3.456664796984375*^9}, {
  3.536606154070375*^9, 3.536606154336*^9}, {3.551717264712678*^9, 
  3.5517172690419254`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"dopplerDist", "[", "v_", "]"}], ":=", 
   RowBox[{"Exp", "[", 
    RowBox[{"-", 
     FractionBox[
      SuperscriptBox["v", "2"], 
      SuperscriptBox["\[Alpha]", "2"]]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{3.551717270225993*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Comb Parameters", "Subtitle",
 FontColor->RGBColor[0, 0.356863, 0]],

Cell["Comb Modes", "Subsubtitle",
 CellChangeTimes->{{3.4566813895*^9, 3.45668139175*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"n1NomA", "=", 
   RowBox[{"Round", "[", 
    FractionBox[
     RowBox[{"\[Nu]GndToInt", "[", 
      RowBox[{"fMinGnd", ",", "fMinInt"}], "]"}], "fRepAvg"], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"n1NomB", "=", 
   RowBox[{"Round", "[", 
    FractionBox[
     RowBox[{"\[Nu]GndToInt", "[", 
      RowBox[{"fMaxGnd", ",", "fMaxInt"}], "]"}], "fRepAvg"], "]"}]}], 
  ";"}]}], "Input",
 CellChangeTimes->{{3.456665147765625*^9, 3.456665238296875*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"n2NomA", "=", 
   RowBox[{"Round", "[", 
    RowBox[{
     FractionBox["1", "fRepAvg"], 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"\[Nu]GndToFin", "[", 
        RowBox[{"fMinGnd", ",", "fMinFin"}], "]"}], "-", 
       RowBox[{"\[Nu]GndToInt", "[", 
        RowBox[{"fMinGnd", ",", "fMinInt"}], "]"}]}], ")"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"n2NomB", "=", 
   RowBox[{"Round", "[", 
    RowBox[{
     FractionBox["1", "fRepAvg"], 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"\[Nu]GndToFin", "[", 
        RowBox[{"fMaxGnd", ",", "fMaxFin"}], "]"}], "-", 
       RowBox[{"\[Nu]GndToInt", "[", 
        RowBox[{"fMaxGnd", ",", "fMaxInt"}], "]"}]}], ")"}]}], "]"}]}], 
  ";"}]}], "Input",
 CellChangeTimes->{{3.456665147765625*^9, 3.4566653128125*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"n1List", "=", 
  RowBox[{"Union", "[", 
   RowBox[{"Join", "[", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"n1NomA", "+", "\[Delta]n"}], ",", 
       RowBox[{"{", 
        RowBox[{"\[Delta]n", ",", 
         RowBox[{"-", "3"}], ",", "3"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"n1NomB", "+", "\[Delta]n"}], ",", 
       RowBox[{"{", 
        RowBox[{"\[Delta]n", ",", 
         RowBox[{"-", "3"}], ",", "3"}], "}"}]}], "]"}]}], "]"}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.456665383359375*^9, 3.45666546196875*^9}, {
  3.456665618125*^9, 3.456665620609375*^9}, {3.5372066701081247`*^9, 
  3.5372066730768747`*^9}, {3.5377187815151253`*^9, 3.5377187956088753`*^9}, {
  3.5377946040932503`*^9, 3.5377946077338753`*^9}, {3.5377949385776253`*^9, 
  3.5377949450151253`*^9}, {3.537893207650926*^9, 3.537893217432143*^9}, {
  3.5385234618671618`*^9, 3.53852346744048*^9}, {3.538583211301021*^9, 
  3.5385832213905983`*^9}, {3.5407344378562293`*^9, 3.5407344440115814`*^9}, {
  3.5409093630926204`*^9, 3.540909373629223*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "1563369", ",", "1563370", ",", "1563371", ",", "1563372", ",", "1563373", 
   ",", "1563374", ",", "1563375", ",", "1565216", ",", "1565217", ",", 
   "1565218", ",", "1565219", ",", "1565220", ",", "1565221", ",", 
   "1565222"}], "}"}]], "Output",
 CellChangeTimes->{3.5521685429609375`*^9, 3.552753835763672*^9, 
  3.552757282147461*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"n2List", "=", 
  RowBox[{"Union", "[", 
   RowBox[{"Join", "[", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"n2NomA", "+", "\[Delta]n"}], ",", 
       RowBox[{"{", 
        RowBox[{"\[Delta]n", ",", 
         RowBox[{"-", "4"}], ",", "4"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"n2NomB", "+", "\[Delta]n"}], ",", 
       RowBox[{"{", 
        RowBox[{"\[Delta]n", ",", 
         RowBox[{"-", "4"}], ",", "4"}], "}"}]}], "]"}]}], "]"}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.456665383359375*^9, 3.45666546196875*^9}, {
  3.456665585125*^9, 3.45666558903125*^9}, {3.456665624328125*^9, 
  3.456665626734375*^9}, {3.5372066753268747`*^9, 3.5372066804518747`*^9}, {
  3.5377187866401253`*^9, 3.5377187978745003`*^9}, {3.5377946088901253`*^9, 
  3.5377946116245003`*^9}, {3.5377949469370003`*^9, 3.5377949552338753`*^9}, {
  3.5378932198813477`*^9, 3.537893225029357*^9}, {3.538523470201638*^9, 
  3.5385234737538414`*^9}, {3.538583224257762*^9, 3.5385832272179317`*^9}, {
  3.540734445712679*^9, 3.540734450881974*^9}, {3.540909375727343*^9, 
  3.5409093833457785`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "1721306", ",", "1721307", ",", "1721308", ",", "1721309", ",", "1721310", 
   ",", "1721311", ",", "1721312", ",", "1721313", ",", "1721314"}], 
  "}"}]], "Output",
 CellChangeTimes->{3.5521685436171875`*^9, 3.552753836354492*^9, 
  3.5527572822109375`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Comb Equation", "Subsubtitle",
 CellChangeTimes->{{3.456681410328125*^9, 3.456681413015625*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"\[Omega]L", "[", 
    RowBox[{"n_", ",", "fR_", ",", "f0_"}], "]"}], ":=", 
   RowBox[{"2", " ", "\[Pi]", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"n", " ", "fR"}], "+", "f0"}], ")"}]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.456678515171875*^9, 3.45667856759375*^9}, {
   3.4566795630625*^9, 3.45667956609375*^9}, 3.5413431977628536`*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Integrate 2-photon Lineshape over Doppler Distribution", "Subtitle",
 FontColor->RGBColor[0, 0.356863, 0]],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"transProb", "[", 
     RowBox[{"\[Omega]1_", ",", "\[Omega]2_", ",", "v_"}], "]"}], "=", 
    RowBox[{
     UnderoverscriptBox["\[Sum]", 
      RowBox[{"fDPrime", "=", "fMinFin"}], "fMaxFin"], 
     RowBox[{
      UnderoverscriptBox["\[Sum]", 
       RowBox[{"f", "=", "fMinGnd"}], "fMaxGnd"], 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          FractionBox[
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{"\[Gamma]Fin", "+", "\[Gamma]T"}], ")"}], "2"], "4"], 
          "/", 
          RowBox[{"(", 
           RowBox[{
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{
               RowBox[{"2", " ", "\[Pi]", " ", 
                RowBox[{"\[Nu]GndToFin", "[", 
                 RowBox[{"f", ",", "fDPrime"}], "]"}]}], "-", 
               RowBox[{"\[Omega]1", 
                RowBox[{"(", 
                 RowBox[{"1", "+", 
                  FractionBox["v", "c"]}], ")"}]}], "-", 
               RowBox[{"\[Omega]2", 
                RowBox[{"(", 
                 RowBox[{"1", "-", 
                  FractionBox["v", "c"]}], ")"}]}]}], ")"}], "2"], "+", 
            FractionBox[
             SuperscriptBox[
              RowBox[{"(", 
               RowBox[{"\[Gamma]Fin", "+", "\[Gamma]T"}], ")"}], "2"], "4"]}],
            ")"}]}], ")"}], 
        FractionBox["1", 
         RowBox[{
          RowBox[{"2", "f"}], "+", "1"}]], 
        RowBox[{
         UnderoverscriptBox["\[Sum]", 
          RowBox[{"m", "=", 
           RowBox[{"-", "f"}]}], "f"], 
         RowBox[{
          UnderoverscriptBox["\[Sum]", 
           RowBox[{"mDPrime", "=", 
            RowBox[{"-", "fDPrime"}]}], "fDPrime"], 
          SuperscriptBox[
           RowBox[{"Abs", "[", 
            RowBox[{
             UnderoverscriptBox["\[Sum]", 
              RowBox[{"fPrime", "=", "fMinInt"}], "fMaxInt"], 
             RowBox[{
              UnderoverscriptBox["\[Sum]", 
               RowBox[{"mPrime", "=", 
                RowBox[{"-", "fPrime"}]}], "fPrime"], "\[IndentingNewLine]", 
              RowBox[{
               RowBox[{"(", 
                RowBox[{
                 FractionBox[
                  RowBox[{"\[Gamma]Int", "+", "\[Gamma]T"}], "2"], "/", 
                 RowBox[{"(", 
                  RowBox[{
                   FractionBox[
                    RowBox[{"\[ImaginaryI]", " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Gamma]Int", "+", "\[Gamma]T"}], ")"}]}], "2"], 
                   "+", 
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"\[Omega]1", 
                    RowBox[{"(", 
                    RowBox[{"1", "+", 
                    FractionBox["v", "c"]}], ")"}]}], "-", " ", 
                    RowBox[{"2", " ", "\[Pi]", " ", 
                    RowBox[{"\[Nu]GndToInt", "[", 
                    RowBox[{"f", ",", "fPrime"}], "]"}]}]}], ")"}]}], ")"}]}],
                 ")"}], 
               SubscriptBox["dGndToInt", 
                RowBox[{"f", ",", "fPrime", ",", "m", ",", "mPrime"}]], 
               SubscriptBox["dIntToFin", 
                RowBox[{
                "fPrime", ",", "fDPrime", ",", "mPrime", ",", 
                 "mDPrime"}]]}]}]}], "]"}], "2"]}]}]}], ")"}]}]}]}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"fnToIntegrate", "[", 
     RowBox[{"fR_", ",", "fOffset_", ",", "v_"}], "]"}], "=", 
    RowBox[{
     UnderoverscriptBox["\[Sum]", 
      RowBox[{"n1", "=", "1"}], 
      RowBox[{"Length", "[", "n1List", "]"}]], 
     RowBox[{
      UnderoverscriptBox["\[Sum]", 
       RowBox[{"n2", "=", "1"}], 
       RowBox[{"Length", "[", "n2List", "]"}]], 
      RowBox[{
       RowBox[{"dopplerDist", "[", "v", "]"}], 
       RowBox[{"transProb", "[", 
        RowBox[{
         RowBox[{"\[Omega]L", "[", 
          RowBox[{
           RowBox[{
           "n1List", "\[LeftDoubleBracket]", "n1", "\[RightDoubleBracket]"}], 
           ",", "fR", ",", "fOffset"}], "]"}], ",", 
         RowBox[{"\[Omega]L", "[", 
          RowBox[{
           RowBox[{
           "n2List", "\[LeftDoubleBracket]", "n2", "\[RightDoubleBracket]"}], 
           ",", "fR", ",", "fOffset"}], "]"}], ",", "v"}], "]"}]}]}]}]}], 
   ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"vIntSpec", "[", 
    RowBox[{
     RowBox[{"fRep_", "?", "NumericQ"}], ",", 
     RowBox[{"fO_", "?", "NumericQ"}]}], "]"}], ":=", 
   RowBox[{"NIntegrate", "[", 
    RowBox[{
     RowBox[{"fnToIntegrate", "[", 
      RowBox[{"fRep", ",", "fO", ",", "v"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"v", ",", 
       RowBox[{"-", "60000"}], ",", "60000"}], "}"}], ",", 
     RowBox[{"MinRecursion", "\[Rule]", "3"}], ",", 
     RowBox[{"MaxRecursion", "\[Rule]", "10"}], ",", 
     RowBox[{"AccuracyGoal", "\[Rule]", "4"}], ",", 
     RowBox[{"PrecisionGoal", "\[Rule]", "2"}]}], "]"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.538523548469115*^9, 3.5385235524603434`*^9}, {
   3.5385832601238136`*^9, 3.538583261230877*^9}, {3.540909398014618*^9, 
   3.540909411994417*^9}, 3.541343202952861*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Calculation", "Subtitle",
 CellChangeTimes->{{3.456681446140625*^9, 3.45668144753125*^9}, {
  3.5379007941050625`*^9, 3.537900798020669*^9}},
 FontColor->RGBColor[0, 0.356863, 0]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"rFrep", "=", "0.0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"tt", "=", 
    RowBox[{"AbsoluteTime", "[", "]"}]}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{
   "distribute", " ", "NIntegrate", " ", "among", " ", "CPU", " ", "cores"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"specList", "=", 
   RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"SetSharedVariable", "[", "specList", "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ParallelDo", "[", 
   RowBox[{
    RowBox[{"AppendTo", "[", 
     RowBox[{"specList", ",", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{
          "fRepList", "\[LeftDoubleBracket]", "kk", "\[RightDoubleBracket]"}],
           "-", "rFrep"}], ")"}], ",", 
        RowBox[{"Chop", "[", 
         RowBox[{"vIntSpec", "[", 
          RowBox[{
           RowBox[{
           "fRepList", "\[LeftDoubleBracket]", "kk", 
            "\[RightDoubleBracket]"}], ",", 
           RowBox[{
           "fOffsetList", "\[LeftDoubleBracket]", "kk", 
            "\[RightDoubleBracket]"}]}], "]"}], "]"}]}], "}"}]}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"kk", ",", 
      RowBox[{"Length", "[", "fRepList", "]"}]}], "}"}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ttt", "=", 
   RowBox[{
    RowBox[{"AbsoluteTime", "[", "]"}], "-", "tt"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"Print", "[", 
  RowBox[{"\"\<Execution time = \>\"", ",", "ttt"}], "]"}]}], "Input",
 CellChangeTimes->{{3.5407346692514644`*^9, 3.5407346703575277`*^9}, {
  3.5413432101228704`*^9, 3.5413432306928997`*^9}, {3.5425568547897897`*^9, 
  3.542556855210992*^9}, {3.5492478888847656`*^9, 3.549247903041992*^9}, {
  3.552168231367849*^9, 3.5521682316408644`*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Execution time = \"\>", "\[InvisibleSpace]", 
   "1759.4179687`10.696914016544744"}],
  SequenceForm["Execution time = ", 1759.4179687`10.696914016544744],
  Editable->False]], "Print",
 CellChangeTimes->{3.552170227366211*^9, 3.5527555276621094`*^9, 
  3.5527590422626953`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"ListPlot", "[", 
  RowBox[{
   RowBox[{"Sort", "[", "specList", "]"}], ",", 
   RowBox[{"PlotRange", "\[Rule]", "All"}], ",", 
   RowBox[{"Joined", "\[Rule]", "True"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.5407062102656565`*^9, 3.5407062108584414`*^9}, 
   3.5407346912477226`*^9, {3.5407454061664057`*^9, 3.54074540677481*^9}, 
   3.540917370755953*^9, {3.540955271524134*^9, 3.5409552734745717`*^9}, {
   3.5409573772453723`*^9, 3.5409573788055224`*^9}, 3.5413432386329107`*^9}],

Cell[BoxData[
 GraphicsBox[
  {Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJwl2Hk8VekDx/GDIqJIUVckpYRENVJq7GlRREWKilRaFEpCuqVpRpQsadNo
HQltPy2mVIQQt4VRSQhZSyltGP3m+61/er1fH/fe55zzPGcb6bXJ2UdaEARL
KUHA//jnU5phFaokfbpn9T0LQRDTfmUfl20LhwUpuPm5klSfKHZ62f96b7sf
Z5eGg8c/eLXvPDsd8fmr4ub17DLwlfAWHY1ydnr++Oze0rvsfeA4B8NHS2Rz
0OlkY495V4xgoS/80jdLvnEuO23f+snR04BdFm5zuanls5SdDh26RdX/ELsc
bD1UI9BWOhedFr1MHlRrBwv94JyGGwvvRrHT5T73O0Mk7PJwRbRZQL3WfXR6
dcPFrWHLYEEBltQ/WuS1kf2nn+5rmvCMvT98ovDYldTleeh0iEHSw2Hz8tEV
4U1THccl1sNiuvfFkcWdGwrQleDpPSeD2s0eoNNaX19unTS7EH0A/F1V99nW
iUXodOWnT6+PbC5GHwiXLE606xxQgk673vQWNT0uRVeGx9W0/lvU8AidrpcM
/yvS9Am6Cjz06oWyh21P0WmV68pmTVP+QR8E3y7ybnZJeY5O3www+n2FzSt0
VbhfyK60C99eo9OlTtqq2aWN6INhz1Nhtk+j3qHTIrVZphXdH9GHwOEiw8Ep
c7+i0yuWi9Sd+nWhq8G6X4ad768Ni+m4ZN1uSQ7+XlCHd6cLj7Z9/oROO34x
TfUb9QF9KMefojb09LxWdHrldpV1mtYc3zD4fMqWWTvW16PTv7udijJPrEUX
wYXjv7l4BMJiOq0hXuv0gmp0Dfjt+BaHiDrsDzEtF2bWbujMPhyeadJ91Emr
Bp0+fser0jyI368Je78ZExViz/1Hp0/cFbHxCyxocT6Ey01I8+b46MWPQ6PX
/tuAPoLr4zeHVpuUJnR6eXeYsc0/2F5Bm/vPTv1+rdx7dLrw4J7OTA/u/5Hw
nwWiwb85fUGnjUbkRbi+7kbXgQv2XS5e+0bK8r9Ou62U2jfbRO4/C6PgudrJ
UkOGD0CnFyVuOOwXMQh9NBz/+XbPzEOD0WntBfoph0JgQRdWriyZ1LgGfy+m
p2oWf5gRiO8TxsBF5RWR2wfJo9Otm+tlnyX3QR8L/yVxb417zfHR6z9c+Lc5
SkDXg/0XPbj636fRaXcjlatV79nHwbZ+hgPHSvh5Om7HgHH+9dLo+rDN9+vr
PJfi98T0i4ze9wln+qIbwDcaohti42ExbWK/3Sa4UgbdkMfLrvbh5If8fnqG
b/+BCS78/fHwwMIURf1FP7D/6QMRJzX+joEFI7hx3N2elrMcP5110iVLfTbH
N4Hn0y9iY7N6jo+OvhzV7pLTD90Y1vFy7Vjjx+NDN+jJDFZ3U0U3gTvVtY7U
6gxDp2uV05WTLLXRJ8LPXp/b5/hyFDqtWV1wL6NXF30SHDspxzmsEBbTlQ9O
1QTvx98Lkzl/kgybE6zwfWL6rX2970VLEfovsF5vlXrGriHotPTC+kVfRnB8
pnDCGZuLbT0q6HTerZtX/LQ4v6bAuRLVPjXG+Hsx7exy/MDbu5xfZnDMu8lL
hH1q6HTkOo9beaFD0afCxi9GpqVGYjxi+nn7qcjU1xzfNDh7Rm7I/NP4ezGd
uvye41dFjFcwh7dYvcvyU+T8pb3dUlbmFg5En87j72ldMcZfCZ2+2fSmXFEJ
FmbwfF1sqDxehceH9tsbc9bWURn9V3jnpW0e4cf4/bTMsdxz96ZgewQLeHSU
4fKJJzTQ6ernCiLfkyPRLXn+kEl4lhI0Bp1u7zCrddUwQLeCBdG2OfWXx6PT
+/VOPHn4CyxY83xgPCi+bqQ+Oh0bIx9z4RGPvw3XX8uZYL8c/J6YrrLriHFo
10S3hZdapFfdrub+pTP2KqlMKeL+tYMt0/anHTXg+OndFTcqotv4+ZlczzNO
56Xlc/7Q9ZGGK+/H6KDbw6F2R2e0rR2NTi+xcGxODMT2CrPgAVlW4xIvjUWn
A24cNbXUYJ8N/9OwMDgxBd8npksz56i81NFCn8P1cGmUv/ZlHn9657GmMJ2R
PP5zef2ft7PyrjXnH118oSSptC/nrwO88Vi7WeIuHj9aFOSn0bGb83Me1/uK
EeoBIZyftJthn6XZ2dw/83l9PZxplxfK7afNfxySOx+B7RUc4Q+PVgdtlddD
p/30v+cuU8XxEpy4/uTnRz2IGodO5y23fr+nltu/gNe7QV/nmphxfdOXkqcp
nNHj+neGjzWKEmqlh6PTNwzfxaZ04XwhuHD/Gp0O8OrD40ufURj+wrkLx1NY
CNekbU5TyBiBTk+75DZ2XBa/fxHHE5KTl3+Cx4/2vbaors8zjm8x16vR9/eZ
1hw/nWvm0zbhV85fV3iTzFp1ryOcn7RcT0ajbTk/7wafqlVwmGDA40s/c+0b
G5uJ7RGWcL2nX0lvGKOOTlsP+rW5MIjHz53reYj1k6YpPP/QWzzb+29Q5vpc
yu23unH4+RWsdzGtVtCztyiefRl8blDqwf7LefzpXtED2Uhvzh8POGqiQn2e
D8+/9Jo7Gbo5Bpz/nrDinQPx59swH8W0p4/qiouynJ/LuT5M60PeSXP90BJj
83rJMGyPsILb1xU1Q/8O5yf9JDXsZlQ6tkdYyevzBoU3Pm8xXjE9y0rvpbE/
uxf8aYJHRK4JP0/rXDZpMjXk+ccb1s8suR+rxuNPf6hscSkYz/Gtgo9YBiUk
NXH+0o4Sac/jV7h+ffj3ucsmVfjx+NPfVudnd2Zyfq/m9bJ/nUHJOh4/+kee
72zZ3dw/a2D1wyerPe9y/dDGQQfvBGZy/Gvh6QbeTmqbFdHp9sOFFn/k8/7F
l9eLfs5PG1uk0ek321uybGV4/V7H6336dNclLux0mfb8laG5sLCe16uVg6t8
NXn9pfclTbXYM0kWfQPcJKdquX8ersdiuv8PqSelQ3A/I2yEL0iMxKLB7PSL
rqAurxLeX/jB6UWhIvvJGI+Ydr5Y0PeOIu4nhE2cz7pXb71Zz/sHWnj81DNP
lvcPmzn/rN7pnzPl/QMdPvD2gpZWft6f+3NZ2L/Tijh+euLxuJ3Vf/H+IQB+
f6nL1LWS+4/+Jctpx5I5nN+B/H7x7OrRnVwftHJX1erhGVw/W+DEbzFORvt5
/0dPnj06tkyOfSv3V35R+PwFXD/00Y7P0wPG8PeD4FWHDEtv+XP/02fD1wYq
3urB9m3j/pucH9LZh88PtJBVFa/SD88DQjAsdh9uclYDzwNiulj+wcQdzW3o
27meUnMKr6c1o9ObPczL1c34fBDC4zU3vzOslvfftJ6VS6diXR16KHyybLHK
1/u8f6cbV+2USpfi/XsYPMrMr8frAJ8faGmHmhVFAbz/3wGPTfBUbTvFTsvb
ev31fRo/H8793V4hOjgHvyemt7m7vopT4v3/Tl6vfWe52odgvGLaplV+b4I5
7//FsI9u/AdrRd7/02OKGmdEevH+fxc8uqM8+9uHTnT6wsMdnbMU+Xy1G167
qbM36e039N0/92fFufJoPp9FwB1PQ4dJ3/qOThdZ5WWq9n5G3wNXeRUqb1vc
gU4PqZm+ZX0Rx/cbzz82c7RlWrA9Ylpq7+DozhLun71wXefHBLWHVei0Ql3S
lN4ReB4Vfuf+l6wty3PH86qYXmZoXdMypAz9D3jENYeSPWp43hXTq1+dX9Mt
wfOwEAkfeqb5IjYZz8ti+tVRT9mpffE8Lezj+fB/nn0fd/P5m558y8FFQQQL
UTze/u5raj7w+Z1OTX5abd/N5/touKW1su/HebCYjtqqaD7pCt8P7Icf+v84
+E6B7wfoi9Wn1P7WxfsE4QCckRXufUma7xdoO8/fSmwc+P4hBr4eEL3nuTHf
T9B7Tq45IXnO9xsHYdkg99F7Cvj+g/7jk9SoMzPZY+FJ2w83P/7M9yu0fbNa
2dW/+X4ljtejVZfraoPZ6ZKPUz85a7LHw0pzrs06+oXvf+gU1/X9tkv4/ieB
56czmjLbz7LT/lFuuuID7IfgCX8WxxsFs9MO18xu/IhkT4Td5aYMt9/KThfH
KU0Y5sR+GLbI3ti1eCY73Xg5/+ZCc/YjcLC6TqCdBTutf2DbRkVT9qPwCUmv
isSNnR55T0lN34T9GH9P6X1poDM7nZ8tJdU+h/043PzdxshzCTs9U9rBdp0r
exK8S/24ZeNGdlp2rOzT4kPsJ3g+0xNdP/lz++mYkNGVcT/3z59wdFhFitTP
9290+Q297pp89mSevwuUOtY3sNNljUEe547cs/g/DioBeg==
    "]]},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  PlotRangeClipping->True]], "Output",
 CellChangeTimes->{3.552170227397461*^9, 3.5527555277685547`*^9, 
  3.5527590427285156`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Export", "[", 
  RowBox[{"oFile", ",", 
   RowBox[{"Sort", "[", "specList", "]"}], ",", "\"\<Table\>\""}], 
  "]"}]], "Input"],

Cell[BoxData["\<\"C:\\\\Users\\\\work\\\\Data\\\\output\\\\K39-D2-4D52.txt\"\>\
"], "Output",
 CellChangeTimes->{3.552170227694336*^9, 3.552175407616211*^9, 
  3.552755528203125*^9, 3.552759043203125*^9}]
}, Open  ]]
}, Open  ]]
},
WindowSize->{1272, 684},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
FrontEndVersion->"8.0 for Microsoft Windows (32-bit) (November 7, 2010)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[545, 20, 486, 13, 72, "Input"],
Cell[CellGroupData[{
Cell[1056, 37, 461, 10, 85, "Input"],
Cell[1520, 49, 377, 10, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[1934, 64, 132, 2, 28, "Subsubtitle"],
Cell[2069, 68, 2015, 33, 31, "Input"],
Cell[4087, 103, 657, 18, 29, "Text"],
Cell[4747, 123, 4468, 92, 332, "Input"],
Cell[9218, 217, 660, 19, 29, "Text"],
Cell[9881, 238, 2672, 75, 314, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[12590, 318, 76, 1, 49, "Subtitle"],
Cell[12669, 321, 37, 0, 28, "Subsubtitle"],
Cell[12709, 323, 191, 6, 31, "Input"],
Cell[CellGroupData[{
Cell[12925, 333, 40, 0, 28, "Subsubtitle"],
Cell[12968, 335, 223, 8, 31, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[13228, 348, 107, 1, 28, "Subsubtitle"],
Cell[13338, 351, 330, 10, 31, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[13717, 367, 75, 1, 49, "Subtitle"],
Cell[13795, 370, 102, 1, 49, "Subtitle"],
Cell[13900, 373, 149, 2, 28, "Subsubtitle"],
Cell[14052, 377, 357, 8, 47, "Input"],
Cell[CellGroupData[{
Cell[14434, 389, 105, 1, 28, "Subsubtitle"],
Cell[14542, 392, 406, 10, 47, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[14985, 407, 101, 1, 28, "Subsubtitle"],
Cell[15089, 410, 219, 6, 31, "Input"],
Cell[15311, 418, 167, 5, 31, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[15527, 429, 106, 1, 49, "Subtitle"],
Cell[15636, 432, 222, 3, 28, "Subsubtitle"],
Cell[15861, 437, 309, 7, 47, "Input"],
Cell[CellGroupData[{
Cell[16195, 448, 199, 3, 28, "Subsubtitle"],
Cell[16397, 453, 285, 7, 47, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[16719, 465, 109, 1, 28, "Subsubtitle"],
Cell[16831, 468, 630, 13, 52, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[17498, 486, 389, 5, 28, "Subsubtitle"],
Cell[17890, 493, 1129, 20, 31, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[19056, 518, 413, 6, 28, "Subsubtitle"],
Cell[19472, 526, 589, 10, 31, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[20110, 542, 356, 10, 52, "Subtitle"],
Cell[20469, 554, 249, 3, 28, "Subsubtitle"],
Cell[20721, 559, 328, 7, 47, "Input"],
Cell[CellGroupData[{
Cell[21074, 570, 199, 3, 28, "Subsubtitle"],
Cell[21276, 575, 463, 9, 47, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[21776, 589, 115, 1, 28, "Subsubtitle"],
Cell[21894, 592, 347, 9, 52, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[22278, 606, 252, 3, 28, "Subsubtitle"],
Cell[22533, 611, 768, 18, 47, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[23338, 634, 315, 4, 28, "Subsubtitle"],
Cell[23656, 640, 1307, 24, 54, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[25000, 669, 350, 5, 28, "Subsubtitle"],
Cell[25353, 676, 1168, 22, 33, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[26558, 703, 371, 5, 28, "Subsubtitle"],
Cell[26932, 710, 1117, 21, 33, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[28098, 737, 538, 12, 52, "Subtitle"],
Cell[28639, 751, 240, 3, 28, "Subsubtitle"],
Cell[28882, 756, 403, 8, 47, "Input"],
Cell[CellGroupData[{
Cell[29310, 768, 199, 3, 28, "Subsubtitle"],
Cell[29512, 773, 656, 12, 47, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[30205, 790, 104, 1, 28, "Subsubtitle"],
Cell[30312, 793, 641, 13, 52, "Input"],
Cell[30956, 808, 827, 20, 75, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[31820, 833, 293, 4, 28, "Subsubtitle"],
Cell[32116, 839, 1022, 20, 47, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[33187, 865, 139, 1, 49, "Subtitle"],
Cell[33329, 868, 175, 2, 28, "Subsubtitle"],
Cell[33507, 872, 442, 15, 31, "Input"],
Cell[33952, 889, 1602, 47, 56, "Input"],
Cell[35557, 938, 1508, 31, 52, "Input"],
Cell[CellGroupData[{
Cell[37090, 973, 194, 4, 52, "Input"],
Cell[37287, 979, 188, 2, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[37512, 986, 57, 0, 28, "Subsubtitle"],
Cell[37572, 988, 1384, 45, 76, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[38993, 1038, 97, 1, 28, "Subsubtitle"],
Cell[39093, 1041, 1254, 33, 99, "Input"],
Cell[40350, 1076, 1184, 32, 99, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[41571, 1113, 115, 1, 28, "Subsubtitle"],
Cell[41689, 1116, 2134, 54, 275, "Input"],
Cell[43826, 1172, 2166, 55, 275, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[46041, 1233, 89, 1, 49, "Subtitle"],
Cell[46133, 1236, 856, 20, 74, "Input"],
Cell[46992, 1258, 352, 9, 53, "Input"],
Cell[47347, 1269, 288, 9, 50, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[47672, 1283, 73, 1, 49, "Subtitle"],
Cell[47748, 1286, 90, 1, 28, "Subsubtitle"],
Cell[47841, 1289, 507, 15, 87, "Input"],
Cell[48351, 1306, 833, 25, 87, "Input"],
Cell[CellGroupData[{
Cell[49209, 1335, 1121, 25, 31, "Input"],
Cell[50333, 1362, 383, 8, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[50753, 1375, 1168, 26, 31, "Input"],
Cell[51924, 1403, 300, 7, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[52261, 1415, 102, 1, 28, "Subsubtitle"],
Cell[52366, 1418, 394, 10, 31, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[52809, 1434, 112, 1, 49, "Subtitle"],
Cell[52924, 1437, 5254, 138, 257, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[58215, 1580, 185, 3, 49, "Subtitle"],
Cell[CellGroupData[{
Cell[58425, 1587, 1856, 52, 152, "Input"],
Cell[60284, 1641, 327, 7, 23, "Print"]
}, Open  ]],
Cell[CellGroupData[{
Cell[60648, 1653, 508, 9, 31, "Input"],
Cell[61159, 1664, 4614, 79, 234, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[65810, 1748, 150, 4, 31, "Input"],
Cell[65963, 1754, 204, 3, 30, "Output"]
}, Open  ]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)
